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Abstract 

The purpose of the present study is to compare the predictions of different models of star 
formation rate (SFR) history in the universe with the upper limit of Super Kamiokande for the 
neutrino background. To this aim we have calculated the expected neutrino density for the most 
popular models of SFR history, Hogg et ah, Glazebrook et ah, Cole et ah, Yuksel et ah, Hernquist 
et al. and Kaplinghat et al. Different from previous studies we have used the ACDM model with 
Da = 0.7. We have assumed that the detector used for the detection the neutrino flux is SuperK 
and also we have assumed that the electron neutrinos produced in the Supernovae oscillate equally 
to the three standard model flavors. By these assumptions all models stay below the upper limit 
of SuperK on the event rate and the detection of the supernova relic neutrino background (SRNB) 
remains undetected. Future neutrino detectors such as KMSNet will be able to detect the SRNB 
and distinguish between the models of the SFR history. 
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I. INTRODUCTION 


Neutrino telescopes are right now the only possibility of exploring the distant universe 
independent of photons. By studying the relic neutrinos stemming from distant supernovae 
we can derive information about the history of the universe and its star formation rate (SFR). 
It is a well established fact that core collapse supernovae like Type II SN emit 99% of their 
energy in the form of neutrinos. These neutrinos, form the so called relic background, 
which is created by Type II supernovae and has been detected in large Neutrino detectors 
such as SuperK, SNOLlab, Icecube and future detectors such as KmSNet. The relic Vf. flux 
depends on the rate of SNu rate as function of redshift which itself is a function of the SFR 
history of the universe. Also the metal enrichment history and the cosmic chemical evolution 
are closely related to SFR and play an important role in the study of the universe. In the 
recent past a wealth of methods for estimating the star formation rate history have been 
applied. For a recent review of the observational methods used in this topic, cf. Madau, 
Wilkinson 2014 [T]. There is a general consensus that the star formation increases from now 
to redshift unity and decreases exponentially towards higher redshifts and earlier times. 

The star formation history of the universe was the subject of a pioneering work by Madau 
et al.96 [2]. Sicne then a large number of studies have been undertaken until the recent past. 
For example Pei and Fall compared models of cosmic chemical evolution with observational 
indicators, the abundance of neutral hydrogen, heavy elements and the dust on damped Lya 
systems and present-day galaxies ESI- Totani et. al EE] used a time-dependent supernova 
rate from a model of galaxy evolution based on the population synthesis method. In a similar 
way, in the work of Bisnovatyi-Kogan and Seidev [7] , galaxy evolution is considered and it is 
assumed that the supernova rate depends redshift proportional to {1 + z)^. In their model, 
the supernova rate is much higher in the early phase of elliptical galaxies, and more than 
half of total supernova explode during the initial 1 Gyr after the formation of galaxies. 
Thereafter the supernovae rate in spiral galaxies becomes dominant and the total number 
of supernova until the present is consistent with the requirements of nucleosynthesis. Also 
Hartman and Woosley [8] compared a model of cosmic chemical evolution with observations 
of Lya systems and faint galaxy surveys. They obtain a a power-law SFR with a redshift 
dependence (SFR oc U^'®|9|). In most models the maximum of SFR takes place at a redshift 
of order unity, where overall rate was about 10 times higher than it is today. 
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In another analysis of the SFR, Hogg considered measurements of radio, infrared and 
ultraviolet broad-band photometric indicators, visible and near-ultraviolet line-emission in¬ 
dicators from redshift unity to the present day uni. Assuming that SFR is proportional to 
(1 + z)^, the best-£t exponent was obtained /? = 2.7 ± 0.7 Cole et. al measured galaxy lumi¬ 
nosity functions in the near infrared from a combined 2MASS-2dFGRS selected galaxy cata¬ 
logue m- On this basis they determined the mass of stars formed until today and obtained 
p* = {a + bz)/[l + {z/cY]hMQyr~^Mpc~^ where(a, 6, c, d) = (0.0166,0.1848,1.9474,2.6316). 
Hernquest and Springel use analytical physical reasons to model SFR.,p* [T2]. They found 
that at early times p* generically rises exponentially as 2 ; decreases, independent of the 
details of the physical model for star formation, but dependent on the normalization and 
shape of the cosmological power spectrum. They conclude that at lower redshifts, the star 
formation rate scales approximately as p* oc H{z)^^^. In this model the peak of SFR de¬ 
pends on the model parameters but half of the stars have formed at redshifts higher than 
2 ; ~ 2.2. Glazebrook and et.al study the overall spectrum of galaxies obtained from the 
red-selected Sloan Digital Sky Survey and compared the results with the blue-selected 2dF 
Galaxy redshift Survay PI- Here they used a double power-law parametrization of the SFR 
with a break at redshift unity: SFR oc (1 -f z)^ for ^ < 1 and oc (1 -1- 2 ;)" for 1 < ^ < 5 
and star formation rate is zero for z > 5.. Finally Yuksel and et al studied SFR in relation 
to the Gamma-ray bursts [H]. The Gamma-ray bursts have the advantage that they can 
be observed at higher redshifts and determine the SFR at z = 4 — 7. The result for SFR 
reported there was that a steep drop exist in the SFR up to at least 2 ; ~ 6.5 

In the present study we have chosen analytic models mentioned above for studying the 
SFR. Here we start out with a model which has a few parameters which can be constrained in 
comparison of such models with observations, which is the relic supernova background in our 
case. This paper is organized as follows. In the next section we calculate the supernova relic 
neutrino density for the different models of SFR. In section 3, we determine the predictions 
of the models for the neutrino event rate at the SuperKamiokande detector. The results are 
presented in section 4, followed by a short discussion in sections. 
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II. THE SUPERNOVAE RELIC NEUTRINO SPECTRUM 


If the supernova rate per unit commoving volume at redshift z is Nsn{Z) and the neutrino 
energy distribution at the source (at energy e) is L^{e), then the expected flux of relic 
neutrinos on Earth is given by [T5] : 


C 


Nsn{z) < LUe') > 


( 1 ) 


(1 + z) a/Oa + (1 + z)^flM 

where e' = (1 + 2 ;)e is the rest frame neutrino energy and the neutrinos are assumed to 
be massless. The spectrum of the neutrinos is parameterized as a Fermi-Dirac distribution 
with zero chemical potential and normalized to the total energy in a particular neutrino 
species {E^) emitted by the supernova, i.e. J L^{e)ede = E^. For each neutrino species Ui 
the energy distribution is given by: 


. ^ 120 e'% , e' 


From the determination of the SFR we extract the SN rate: 


( 2 ) 


Nsn{z) oc p{z) (3) 

where p{z) is the star formation rate. Averaging over a Salpeter Initial Mass Function (IMF) 
for M > 8Mq, the supernova rate is Nsn{z) = while the star formation rate is 

measured in solar masses. Following iMinniiis], we parameterize the SFR as: 


p{z) oc (1 + z)^ (4) 

where the constant A has a different value for z < 1 and 1 < ^ < 2. We also assume that 
the behavior of A at 1 < z < 2, continues to higher redshifts. For z < 1 Hogg [10] has 
compiled measurements of the UV and ffa luminosity density and obtained the 68% C.L 
limits of A = 2.7 ± 0.7 for ACDM model and A = 3.3 ± 0.8 for = 1 model. From the 
optical spectrographic measurements in the Sloan Digital Sky Survey (SDSS) limits on A 
were found to be 2 — 3 for < 1 and 0 — 1 for 2 ; > 1, na. With these assumptions the rate 
of relic neutrinos on earth will be: 

T C 120 <E^> 0.013 f (1 + z)^d^ _ 1 

Ho 7 tt^ <T^ Mq J (1 + 2:)v/DrTTIT^)%A()^exp(|^) + 1 
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( 6 ) 


As before e' = (1 + z)e. By setting x = 1 + z , we can rewrite the integral: 

^ C 120 <E^> 0.013 2 r , 1 

Ho7n^ < Mq a/IIa + x^^Im exp(^) + 1 

The data from S'A^1987A gave = 8 x 10^‘^ergs and = 4.8Mev, while the model 
introduced by Woesley et aim] predicts that for a 25M© supernova progenitor: E^, = 
11 X lO^^er^i, Tj, = 5.3Mev. As an approximate solution for (|^ we hnd: 


C 120 <E^> 0.013 e 
^ >3 Mq Qm 


[1.8^/nZ+^ - O.Sy/f^A + 9flM] 


(7) 


But the SFR is not always simple and it has a more complex shape in the general case. Cole 


et.al obtained from high-z galaxies and gamma-ray bursts data,assuming a Salpeter initial 
mass function HH 


p{z) = {a + bz)h/[l + [z/cY] 


( 8 ) 


with a = 0.0389 ,b = 0.0545 , c = 2.973 and d = 3.655 . Hermquist and Springel found 

m 

pY) = PoxV[l + «(X - 1)^ exp(/3x^/^)] (9) 

where y = [H{z) / , po = 0.030 ,a = 0.323 and /? = 0.051. In the Yuksel et.al model, 
the SFR is [E] 

Mz) = P„l(l + zr + {(1 + + {(1 + ^j/cnvyii + a(x - 1)’ exp(/3,y-'-)| (10) 

where R = (1 -|- C = (1 -|- while po = 0.0285, a = 1.6, 

fd = —1.2, 7 = —5.7 in which zi = 1.7, Z 2 = 5.0 and p = —1.62. Therefore the relation Q 
becomes 

- ^ <Ey> 0.013 2 /■“ , xp(x) _ 1 

Ro77r4<Tj,>4 Mq ^ + x'^VLm exp(^) -h 1 

Neutrino density function for each model has been plotted in hgurel. The results presented 
here in section .4 are based on numerical calculations of ([^ and ( [IT| ) in the intervals z = 0 
to 2 ; = cx) or X = 1 to X = cx). 


III. EVENT RATE AT SUPER KAMIOKANDE DETECTOR 

State of the art neutrino telescopes are water Cerenkov detectors, where two general 
types can be distinguished. The dominant reaction in a light water Cerenkov detector such 
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FIG. 1: Neutrino density predictions at different redshifts for each model as explained section 2. 

as SuperK is PeV —> ne'^ with the cross section ^^(e) two order of magnitude larger than 
that of the scattering reaction Vf.e —)■ in a heavy water Cerenkov detector such as SNO 
is vD —)■ nne'^, where the cross section of these reactions are denoted (jj. z is a positron at 
SuperK type of detector and Deuterium at the SNO type. For simplicity at the calculation 
of the event rate R, the efficiency of the detectors are assumed to be 100% in the observable 
energy window. SuperK does not detect supernovae relic neutrinos at all energies. Below 10 
Mev, the Ue flux is dominated by nuclear reactors and other sources of Pg arriving on Earth 
are not distinguishable. Between 10 Mev and 19 Mev the Neutrino background is due to 
solar neutrinos and due to Cosmic muons in detectors. Between 19 Mev and nearly 20.3Mev 
the background is due to atmospheric neutrinos and the flux of neutrinos at energies greater 
than 36 Mev rapidly falls exponentially, so that the observable flux for neutrinos in SuperK 
detector is from 20.3 to 36.3 Mev where e = Eg + 1.3Men and is the energy of positron 
and 1.3 Mev the neutron-positron mass difference. Therefore the differential event rate in 
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the internal de at SuperK is 

RsuperK = X 1033)(9.52 X 


piz) 




= 5(1.51 X 1033)(9_52 X 1052 ) <^^>J ^) r xp{x)dx 

< Ty Mq Ji y/Q]^x^ + Ha 


X 


X 


.36.3 


l.Syde 


'20.3 




^36.3 g2(g _ 13)2^^ 


'20.3 


exp( 


5f>) + l 

( 12 ) 


where B = {^)CHq ^ = lOSGhgQ^ , we use cTp(e) = 9.52 x lO^^^ePeC'T^^ [IE] ,< E,y >= 
11 X lO^^ergs and < Tj, >= 5.3Mev. So the differential event rate in the interval de is 
NpapJu{e)de and the predicted event rate at the detector are given by: 


R = BN. 


< E^> ,0.013, 


xf>{x) dx 


e^cri(e) 


< >4 ' Mr, 


(13) 


O J V^A + J ®^P(<f^) + 1 
The de delineate the energy window between 20.3 to 36.3 and Np is the number of free 

protons in SuperK detector with a sensitive water mass of 22.5 Ktons and Np = 1.51 x 10^^ 
. The SN relic event rate at SuperK in AC DM model can be written as 

events 


5 = 0.063(y^)( 


Mo (5.) 




(14) 


{Mz) ^ 10^3g7-^s ^ ^ Mev ^ 22.5kton — yaer 
We have set h^o = 1 and also the average metal yield per supernova have taken to be IMq 
in the interest of obtaining an upper bound to the event rate, while the first number in the 
right side of the above expression with Qm = 1 is equal to 0.066. 

More recent estimates by Totani et al [5], using the population synthesis method to 
model the evolution of star formation in galaxies, obtained a prediction for the flux of SRN 
at superK (in the energy interval from 15 to 40 Mev) of 1.2yr~^ and the ’’most optimistic” 
prediction for their model was an event rate of 4.7j/r“^. Malaney [19] used the Pei and Fall 
result [3] in order to parametrize the evolution of the Cosmic gas density rate and integrated 
over all energies and found a total SRN flux of 2.0 — 5.4cm“^s“^. Hartmann and Woosley 
[8] used an SN rate proportional to (1 + z)^ and their best estimation was ~ 0.2cm“^sec“^ 
. Kaplinghat et al. [T3] used the assumption that the supernova rate tracks the metal 
enrichment rate can be written as: Nsn{z) = where < M^ > is the average yield 

of ” metal” {z > 6) per supernova and pz is the metal enrichment rate per unit comoving 
volume 22.bKton — year and the neutrino flux at superK is 1.6cm“^sec“^ or event rate is 
R < 4events for 19MeV < Eg < 35MeV and over all energies the event rate at SuperK is 
54 cm~‘^sec~^. 
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Star Formation Rate 

The neutrino flux (cm ^sec 

Event Rate( R < ) 

Kaplinghat’s model 

1.54 

3 

Hogg’s model 

108.57 

212 

Glazebrook’s model 

111.43 

217 

Cole’s model 

3.02 

6 

Yuksel’s model 

80.32 

156 

Hennguiest’s model 

70.13 

136 


TABLE 1: The neutrino flux and upper bound for SRN at superK in ACDM cosmology 

IV. RESULTS 

Here to obtain the most optimistic of SRN event rate at SnperK we consider nentrino 
oscillation as a mechanism for maximizing the SRN flnx. We have assnmed that AIm = 0.3 
and = 0.7. In Kaplinghat et al model with Ha, our estimates gives the neutrino flux 
at superK is 1.54cm“^sec“^and the event rate is i? < 3 events. For Hogg’s model wich 
has compiled the UV and Ha luminosity density and obtained the SFR as (1 + ^)2.7±o.7 
flux at detector is 108.57 cm~‘^s~^. For the SFR obtained from the optical spectrographic 
measurements in the Sloan Digital Sky Survey (SDSS) the maximum flux of SRN is 111.43 
and minimum flux is 91.1 cm~‘^s~^. With the Cole model the flux is 3.02 
and Yuksel model the flux is 80.32 and then for Hernquest model the flux at detector 

will be 70.13 So then their upper bound from models of Fall, Hogg, Glazebrook, 

Cole, Yuksel and Hernquest are 3 event, 212 event, 217 event, 6 event, 156 event and 136 
event respectively and as summarized in the table [T| 


V. DISCUSSION 

In the present work we have studied the predictions of the different models of SFR for the 
number of expected relic neutrinos in large neutrino detectors such as SuperK. These relic 
neutrinos have been produced at every supernova explosion throughout the history of the 
universe, which transfer almost all their energy to neutrinos. Here we have calculated the 
neutrino production for different SFR models and have found that the models of Glazebrook 






[13] and Hogg |T0] yield an estimate for the nnmber of relic nentrinos closer to the lower 
limit presently given by SnperK. As the relic nentrino backgronnd is an independent way 
of constraining the SFR of the nniverse, this method can be regarded as complimentary 
to the more standard methods. As the constrnction of more modern nentrino telescopes 
is underway, the predictions of the SFR models for the relic neutrinos can be soon tested 
with much more and better data of these telescopes. Regarding this perspective the present 
study can be expected to be extended soon and yield more fruitful results in the light of 
new data. 
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